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Abstract 

We present a new method for studying high-py dilepton events (e^e^, n^fi^, e^/U^) and simulta- 
neously extracting the production cross sections of pp — tt, pp — W^W^, and pp ^ ^ t^t^ 
at a center-of-mass energy of ^/s = 1.96 TeV. We perform a likelihood fit to the dilepton data in a 
parameter space defined by the missing transverse energy and the number of jets in the event. Our 
results, which use 360 pb^^ of data recorded with the CDF II detector at the Fermilab Tevatron 
Collider, are a{ti) = 8.5l|^pb, a{W+W-) = 16.3ltjpb, and a{Z^ t+t') = 29ll^g pb. 

PACS numbers: 14.70.-e, 13.85.Qk, 13.85.Ni 
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There are relatively few standard model (SM) processes that contribute significantly to 
a final-state containing a pair of highly energetic charged leptons. The processes that can 
contribute to these "dilepton" events include top-quark pair (tt) production, VT-boson pair 
production, and Drell-Yan processes. The distinctiveness of this final-state offers 
unique tests of the SM and an intriguing potential for revealing new physics. For instance, 
in Run I (1992-1996) the Collider Detector at Fermilab (CDF) observed several ti candidate 
events in the dilepton decay mode [1] with unusual characteristics, and it was suggested 
that the kinematics of these events could be better described by the cascade decays of heavy 

n 

supersymmetric quarks [2]. Furthermore, the top-quark's extraordinarily large mass might 
be an indication of a close connection with the mechanism of mass generation itself, as for 
example, in the model of topcolor assisted technicolor j^, which predicts new resonances 
decaying to tt. In addition, a fourth generation fermion family ^ could enhance the gluon- 
gluon fusion Higgs production cross section by an order of magnitude, which, for a heavy 
Higgs boson, would lead to an increase in the number of W boson pairs observed |^. These 
examples could all produce an excess of dilepton events in our data, which, depending on 
the topology of the new physics, would affect the cross section measurements of the main 
SM processes to different extents. This provides the motivation behind the present study of 
highly energetic dilepton events. ^ 

The main SM processes with a high-p^ IS] dilepton final-state can be identified based on 
their distinct event characteristics. For example, ti and W^W~ production with decays in 
the e/i dilepton channel, ti —>■ W^hW^h e^fi^uubb and W~^W~ —>■ e^fi^uu, and the 
di-tau decays of bosons, Z° — r+r^ e'^fi^ui/uu, can be distinguished from each other 
by the number of jets, Nj, and the missing transverse energy ^t, in the event. Both 
ti and W~^W~ events typically have large from the final-state undetected neutrinos; 
however, due to the two final state 6-quarks, ti has a greater number of jets than W^W~ . 
Conversely, —>■ t^t~ events have small J^t (due to the neutrinos being of lower energy 
than in the ti and W^W~ processes, and typically traveling in opposite directions), and 
most often no jets. For both W^W~ and Z^ t^t~ events, jets can arise only through 
higher-order processes that include initial-state gluon radiation. 

In contrast to analyses dedicated to measuring a single SM process, the analysis presented 
here adopts a more global strategy by considering all events with a high-py electron and 
muon, and making no further selection requirements. We then exploit the different J^t and 
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Nj characteristics to simultaneously extract the production cross sections of the three main 
processes described. This is done by fitting the efi data in a two-dimensional (2-D) — Nj 
parameter space to template distributions of ti, W^W~, and t^t~ events. Less 

significant processes are also taken into account. 

We also consider the ee and /i/i final-states, in addition to e/i, but in these cases we 
have the added complication of a large Drell-Yan {Z^/'-f* e^e~ or yU^/x^) contribution, 
necessitating a different treatment for these channels. This involves reducing the Drell- 
Yan ee and /i/i contributions by requiring events to have significant Ifix in those channels. 
Without this requirement the ti and W^W~ contributions to the ee and /i/i events would be 
overwhelmed, rendering these final-states unusable. Our ti and W^W~ results use all three 
dilepton final-states. For extracting the Z^ t^t^ cross section we use only e/i events 
as the removal of Drell-Yan ee and /z/i events also significantly reduces the Z^ t^t~ 
contribution in the ee and /i/i channels. 

Since this method makes minimal requirements on events after requiring two leptons, it 
utilizes the full statistical power of the data for given lepton definitions. In the present anal- 
ysis we have chosen to use very tight lepton identification requirements to demonstrate the 
method, and establish a foundation for future measurements with more data. This gives us 
greater control of the background processes, in particular those involving jets being misiden- 
tified as, or containing, leptons. However, it also reduces the potential statistical gain from 
the method. Even so, the results are comparable in precision to the analyses dedicated to 
the individual cross section measurements, and that use looser dilepton definitions. Further- 
more, by looking at all processes simultaneously in the same generic dilepton sample, this 
method tests the SM consistency in a way an analysis dedicated to a single cross section 
measurement does not. An additional benefit of minimizing the requirements on the event 
after two high-p^ leptons are selected is the possibility for sensitivity to new physics which 
might fall into the If^T — Nj parameter space that we study. 

The data sample was collected with the CDF II detector between 2002 and 2004, and 
corresponds to an integrated luminosity of approximately 360 pb~^. The results presentee 
here complement the individual cross section measurements from CDF for each process [a, 
|9|, llOj, which make additional requirements to reduce, and assume the SM cross sections for 
all the processes in the dilepton sample other than the one being measured. 

The CDF II detector has a general-purpose design The components relevant to 
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this analysis are briefly described here. A tracking system inside a 1.4 T superconduct- 
ing solenoidal magnet is composed of silicon detectors for high-precision track measure- 
ments, and an open-cell drift chamber surrounding the silicon system. Electromagnetic and 
hadronic calorimeters surround the tracking system and solenoid and are used to measure 
the energy of interacting particles. In this analysis information is used from both the central 
(covering the pseudo-rapidity range |?7| < 1.1) and end-plug detectors (of which we use the 
range 1.2 < \ri\ < 2.0) to identify electron candidates. The missing transverse energy calcu- 
lation uses calorimeter towers with \ri\ < 3.5. We identify jets using clusters of calorimeter 



towers above an energy threshold of 3 GeV, with fixed cone radius AR = y Acj) + Ar] = 0.4. 
The jet transverse energy is corrected for the calorimeter response and multiple interac- 



tions 
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13| . We require that the corrected jets have Et > 15 GeV and |?7| < 2.5. A set of 



drift chambers located outside the central hadronic calorimeters are used to detect muons 
in the region \ri\ < 0.6. Additional drift chambers and scintillation counters detect muons 
in the region 0.6 < |?7| < 1.0. 

A three-level trigger system is used to select events. The triggers employed to collect 
events for this analysis are |l4] an inclusive central electron (|?7| < 1.1) trigger requiring an 
electron with Et > 18 GeV, and an inclusive central muon (|?7| < 1.0) trigger requiring a 
muon with pt > 18 GeV/c. 

Events selected for the analysis contain two opposite sign leptons (electrons or muons) 
consistent with originating from the same vertex, and with Et > 20 GeV for electrons and 
Pt > 20 GeV/c for muons. Muons from cosmic rays and electrons from photon conversions 
are removed, as described in Ref. llm. Both leptons are required to be isolated in the 

n n 

calorimeter [15[ and the tracking chamber |16|], in order to reduce the probability that a jet 
is misidentified as a lepton, or that a selected lepton comes from the semileptonic decay of 
6 or c hadrons. 

We consider three classes of dilepton events: e/i, ee, and fifi. For the e/i channel we make 
no further requirements on the event after the selection of two leptons as described above, 
and simply count the number of jets, Nj, and measure the J^t in the event. 

The ee and fifi channels have a large Drell-Yan contribution which we reduce significantly 
by applying a requirement on the missing transverse energy significance defined by: 

vJ2 Et 
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where J2 Et is the sum of transverse energies over all calorimeter towers, corrected to include 
the pt of the muons. We require > 2.5 GeV^^^ for all ee and /x/x events . 

After the above dilepton and fi^^ requirements, the dominant SM contributions which 
we consider as our signal processes, are tt — > W~^bW~^b i^ubi'ub, W^W~ —>■ i^ui'u, 
and T^T^ — s> t^viVri'vtVri where I indicates an electron or muon. These processes 

are separated in our chosen If^j- — A'^- parameter space as previously explained. The predicted 
distributions of these signal processes in the ~ parameter space for the e\i channel 
are shown in Fig. [1], and were obtained from Monte Carlo simulations. 

Contributions to the dilepton sample which we consider as background processes, include 
Drell-Yan (ee, /x/i), WZ^ ZZ^ + 7, and W + jets. Note that Drell-Yan events could 
be reconstructed as an e/x final-state when a muon in the forward region radiates a photon 
that is misidentified as an electron. The predicted combined background distributions in the 
IpT — Nj parameter space is also shown in Fig. [TJ We fit the data e/i, ee, and /iyU distributions 
to the expected signal shapes, letting each of their normalizations float in the fit. The e/x 
data distribution is shown in Fig. [21 of which about 60% of the events are expected to be 
from Z^ T^T~ and concentrated in the low-^y, zero-jet region of the ~ Nj parameter 
space. 

To determine the acceptance of our selection criteria for the ti, W^W~, and Z^ t^t~ 

n 

processes we use the pythia p/7| Monte Carlo program, followed by a full simulation of the 
CDF II detector which is based on the geant simulation program The acceptances 
for each of these processes are shown in Table [B The Z^ t^t~ acceptances are defined 
as the fraction of Z^/'y* t^t~ events generated in the di-tau mass range of 66 < M^-t < 
116 GeV/c^ that pass our dilepton selection criteria. The contribution from 7* is only about 
0.3%. For the ee and fifi channels the additional requirement removes 32% of tt events 
and 36% of W^W~ events, while also removing about 99.8% of Drell-Yan events. 

The systematic uncertainties on the ti, W~^W~, and Z^ —>■ t^t'^ acceptances, summa- 
rized in Table Hlj result from uncertainties in the jet energy scale (JES) [l3|, the modeling 
of the initial-state radiation (ISR) by pythia [and for the case of ti also final-state radia- 
tion (FSR)], the uncertainty on the parton distribution functions [19], the modeling of , 
and uncertainties in the lepton trigger and identification efficiencies. The f^^^ systematic 
uncertainty does not apply to the Z^ —>■ t^t^ cross section measurement, as we fit for this 
only in the e/i channel. Moreover, the e/i channel acceptance does not suffer from a jet 
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FIG. 1: The I^t — distributions for tt, W'^W~ , t^t^ in the e/x channel. Also shown is 

the combined background distribution for the e/i channel. All the distributions are normalized to 
an arbitrary equal volume. The numbers of events expected in 360 pb~^ for each source are given 
in Table lllli The highest and Nj bins include any overflow events. 

energy scale uncertainty, as we include efi events with all jet multiplicities, whereas for the 
ee and fifi channels this systematic uncertainty enters through the requirement. A 6% 
uncertainty on the integrated luminosity is applied to the expected number of events for all 



processes \22 \. 



The Drell-Yan and diboson {WZ, ZZ) backgrounds are determined using pythia Monte 
Carlo, followed by the detector simulation. We normalize the total number of events for these 
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e|a, data (103 events) 




FIG. 2: The — distribution for data in the efi channeL The log scale is used to make the 
low count bins more visible. 



TABLE I: Summary of acceptances for tt, W^W', and t'^t~ events, where the quoted 

errors include the systematic uncertainties from Table [Til Values include SM branching fractions 
to the dilepton final-state. The tt events were simulated with a top quark mass of 178 GeV/c^. 



ee 



tt 

w+w- 

Z ^ TT 



(0.399 ± 0.029)% 
(0.294 ±0.018)% 
(0.0458 ± 0.0032)% 



(0.144 ±0.019)% 
(0.111 ±0.008)% 
(0.0008 ± 0.0001)% 



(0.136 ±0.015)% 
(0.092 ± 0.006)% 
(0.0005 ± 0.0001)% 



processes to theoretical cross section predictions 23|]. To estimate the W + '-f background we 
use a matrix element generator j24] and use pythia for the initial-state QCD radiation and 
hadronization. The background from W + jets, where a jet or track is misidentified as an 
electron or muon, is determined from the data. We first calculate the probability that a jet 
with a large fraction of its energy deposited in the electromagnetic calorimeter is misidentified 
as an electron, and the probability that a minimum ionizing track is misidentified as a muon. 
These probabilities are termed fake rates. The fake rate for each lepton type is calculated 
using an average of four inclusive jet samples (triggered with at least one jet with Et > 20, 
50, 70, or 100 GeV respectively). We remove sources of real leptons from Z decays using an 
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TABLE II: Summary of systematic uncertainties on the acceptance for each "signal" process. See 



text for further details. 


Source 


ZL[^ee ) 


++( C 1 l\ 

zzyefi ) 


++(11 1 1\ 


W+ (cc\ 

vv vv (^ee ) 


vv Vv vv vv Z/ — 




JES 


5% 




6% 


1% 


1% 




ISR 


8% 


4% 


6% 


5% 


5% 5% 


5% 


FSR 


7% 


3% 


5% 








Other 


6% 


5% 


5% 


5% 


4% 4% 


5% 


Total 


13% 


7% 


11% 


7% 


6% 7% 


7% 



invariant mass cut, and from W decays using a Monte Carlo estimate of the contamination, 
and parametrize the fake rates as a function of jet transverse energy for electrons, or track 
transverse momentum for muons. The background is determined by weighting the jets from 
a data sample of {W iu) + jets events by the fake rates. As a result of very low statistics 
in the data for the calculation of the fake rates, and large uncertainties in other aspects 
of the calculation, we assume a 100% total uncertainty on our final W + jets background 
estimates. 

A summary of all expected contributions for each dilepton channel is given in Table IIIIl 
together with observed numbers of events. 

We extract the tt, W~^W~, and t^t^ cross sections simultaneously for the e^x 

channel by maximizing the binned likelihood function: 

L{fitt,fj'ww,fj'TT) = n ~ — j — ^ n^i^j^^j) (1) 

where the index i runs over all bins in the two-dimensional Ifix — Nj parameter space, 
and j runs over all variables Xj, which are parameters in the likelihood function constrained 
by a Gaussian of width cij, the estimated uncertainty on Xj. These variables consist of the 
expected number of events for the background processes (given in Table [In|). the acceptances 
for the signal processes (given in Table [T]), and the integrated luminosity. 

The parameters t-i-ww^ and the expected total numbers of tt, W^W , and 

—>■ T^T^ events, respectively. The expected distributions of these processes in the J^t—Nj 
parameter space determine the probabilities ptt,i, Pww,i, and pT-T,i that a ti, W'^W~, or 
Z^ — s> r+r^ event, respectively, will appear in the i-th bin. Therefore, the total expected 
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TABLE III: The numbers of SM predicted events, and the numbers observed, in 360 pb~^ of data. 
For the e/i channel, where the only requirement is two high-pj- leptons, the first three processes 
are considered signals for which cross sections are measured. For the ee and /i/i channels, where 
an additional ^f^^ requirement is made, only the first two processes are regarded as signals. To 
calculate the expected number of events from our signal processes we used the cross section central 
values, cr(tt) = 6.1 pb, a{W^W') = 12.4 pb, and o-(Z° t+t~) = 251 pb. Uncertainties on the 
theoretical cross sections are not included. 





e^u 


ee 






tt 


10.0 ±0.7 


3.6 ±0.5 


3.4 ±0.4 


17.0 ± 1.6 


w+w- 


13.8 ±0.8 


5.2 ±0.4 


4.3 ±0.3 


23.3 ± 1.5 


^ r+r- 


57.8 ±4 


1.1 ±0.2 


0.6 ±0.1 


59.5 ±4.3 


DY ee 





15.4 ±3.2 





15.4 ±3.2 


DY fifi 


9.3 ±0.8 





11.6 ±2.4 


20.8 ±3.2 


WZ 


0.70 ± 0.06 


1.26 ±0.09 


1.11 ±0.08 


3.07 ±0.23 


ZZ 


0.07 ±0.01 


0.47 ±0.03 


0.42 ± 0.03 


0.96 ± 0.07 


W-f 


1.2 ±0.5 


1.8 ±0.7 





3.0 ± 1.2 


W + jets 


3.0 ± 3.0 


2.1 ±2.1 


1.6 ± 1.6 


6.8 ±6.8 


Total SM 


96 ±5 


31 ±4 


23 ±3 


150 ± 12 


Data 


103 


24 


29 


156 



number of events in the z-th bin of the Ifir — Nj parameter space is: 

f^i = {Ptt,i X fitt) + {PWW,i X I^WW) + {PTT,i X firr) + riothev (2) 

where Uother is the expected number of events in the z-th bin from all background processes, 
and is fixed in the likelihood fit within its Gaussian constraint. The product over i in the 
likelihood function is the product of Poisson probabilities for each bin, in which the expected 
number of events is /ij, and rii is the corresponding number of events observed in the data. 
We use 10 GeV wide bins for < 60 GeV, and 20 GeV bins for 60 < I^t < 120 GeV. Note 
that in Figs. [1] and [2] we show the — shapes using a uniform 10 GeV binning in 
For each of the signal processes, fik = o'k^k'^, where k = tt, W~^W~, or Z° — r+r^. 
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and (jfc is the production cross section for the process k, which is a free parameter in the 
hkehhood fit. The acceptances, e^, and integrated luminosity of the data sample, C, are fixed 
within a Gaussian constraint as mentioned above. By maximizing the likelihood function 
we extract the production cross sections for ti, W~^W~ , and t^t~ . 

We perform a similar likelihood fit to the full ee + e/i + /i/i data to extract tt and W^W~ 
cross sections, in which we consider t^t^ as a fixed background included in the 

"Mother term, since in the ee and /i/i channels it has been significantly reduced by the fi^^ 
requirement. We perform the full fit using a likelihood function that is the product of the 
individual likelihood functions for each channel. We also make the following assumptions 
about correlations in the full fit: within a given channel (ee, e/i, or /i/i) we assume the 
signal acceptances are 100% correlated because they are driven by the lepton identification 
efficiencies, and between channels we assume no correlations in the acceptances. The latter 
is not actually true because, for example, the ee and e/i channels have correlations in the 
lepton identification uncertainties due to overlap in lepton types, and the ee and fifi channels 
share some JES systematic uncertainties. However, by varying these correlations between 
0% and 100% we see a negligible effect on the extracted cross sections, so we use 0% for 
simplicity. 

For both the e/t-only and ee + e/t + /t/t scenarios we perform two sets of fits. Our 
main results are obtained by letting all the signal cross sections fioat, with the exception of 
Z^ T^T^ in the ee and /i/t channels, to simultaneously extract the signal cross sections 
from the fit. These results are summarized in Table llVt with the systematic uncertainties 
included in these results being discussed below. As a result of the good separation of our 
signal processes in the fix — Nj parameter space, we observe very little correlation between 
these cross section measurements. In the e/i channel fit where all three cross sections fioat, 
these correlations are about —0.06 between a{tt) and a{W~^W~), —0.05 between a{tt) and 
o-(Z° r+r-), and, -0.19 between a{W+W-) and a{Z^ r+r"). The ti production 
cross section measurement is relatively insensitive to the top mass used for the tt acceptance 
in generating the — Nj template shape. We generated templates using a top mass of 

! 

178 GeV/c^ [25]. Using simulated experiments we observe a 1% variation in a{tt) if the top 
mass is varied between 165 GeV/c^ and 178 GeV/c^; therefore, we neglect any effects due to 
the uncertainty on the top mass. 

Fits to the data are also performed by fixing all but one of the signal processes to their 
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SM expected values. We consider the results extracted from these fits as cross-checks, 
as they use the added constraint of assuming SM production for all processes other than 
the one being measured, in a similar fashion to the more standard counting experiment 
results. Significant differences between these cross-checks and our main results could be 
an indication that the SM assumptions being used are incorrect. We summarize the results 
from these fits in Table IVl When fixed in a particular fit, we use the following SM theoretical 



predictions: a(tt) = 6.1 ±0.9 pb 



251.3 ±5.0 pb 28|. 



26|, a{W+W~) = 12.4±0.8 pb [271], and a{Z^ r+r") 



The uncertainties on the measured cross sections include a component coming from the 
fit (including statistical, acceptance systematic, and integrated luminosity) and a second one 
due to changes in the — Nj distributions caused by the systematic sources mentioned in 



Table [m The fit program used was MINUIT 29|], which minimizes — ln(L), and determines 



the fit errors from the —2 ln(L) values. To evaluate systematic changes in the shape of 
— Nj distributions, we use simulated experiments with one of the signal or background 
distributions from a Monte Carlo simulation with a particular systematic effect applied. The 
shape systematic uncertainties are summarized in Table IVI[ 

TABLE IV: Cross section measurements from a global fit of 360 pb~^ of high-pT dilepton data. The 
first uncertainty is the error returned by the likelihood fit, which includes statistical, acceptance 
systematic, and luminosity uncertainties. The magnitudes of the latter two uncertainties are given 
in Table lU and in the text. The second is the systematic uncertainty in the template shapes. 

Process efi ee + j^ifi + 

^ 9.3l|^ pb 8.5;|^ pb 

a{W+W-) 11.410 ioi pb 16.3;ii pb 



In summary, we present a new method to study globally the production of events with 
final-states including two high-pj' leptons. We measure simultaneously the production cross 
sections for tt, W~^W~, and t^t^ and obtain the following results: a{tt) = 8.5II2 

pb, a{W+W~) = 16.3tli pb, and a{Z^ r+r") = 29lt^|] pb. They are in good agreement 
with SM theoretical predictions. In addition to the potential this analysis technique has for 
precision cross section measurements in the dilepton channel with more data, it could also 
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TABLE V: Cross section measurements from a fit to the data with aU but one signal process fixed 
to its SM value. See text for further details. The uncertainties have the same meaning as in 

Table HVl 

Process efi ee + n/i + 

a{tt) {W+W-, Z° T+T- fixed) 
(j{W+W') {tt, Z° T+T- fixed) 
a{Z tt) {tt, W+W- fixed) 



q 0+3.1 +0.7 u 
^•"5-2.6 -0.2 

12.310 Pb 



i6.iii°i°:« pb 



TABLE VL Summary of the shape systematic uncertainties for the e^u and full(ee + e/i + /i^u) fit. 
Source tt(e/x) W+W-{efi) Z°^T+T-{en) (full) Ty+VF^ (full) 

JES tl% 1:1% ~l% ~l% 

ISR tt% ±1% ±1% 11% ±1% 

FSR ±1% — — ±1% — 

Total tl%o 11% tl%o 



be promising for model independent searches for new physics. 
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